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Abstract. A robust measurement of the clustering amplitude of the sub-mm population of starburst galaxies 
requires large-area surveys (S> 1 deg 2 ). The largest-format arrays subtend only 10 arcmin 2 on the sky and hence 
scan-mapping is a necessary observing mode. Providing realistic representations of the extragalactic sky and 
atmosphere, as the input to a detailed simulator of the telescope and instrument performance, allows important 
decisions to be made about the design of large-area fully-sampled surveys and observing strategies. In this paper 
we present preliminary simulations that include detector noise, time-constants and array geometry, telescope 
pointing errors, scan speeds and scanning angles, sky noise and sky rotation. 



GENERATING REALISTIC SYNTHETIC TIME-SERIES 

A mapping simulator has been developed to generate synthetic bolometer time-series data and associated astrometric 
information which are then run through a realistic reduction pipeline (http://www.inaoep.mx/~echapin/scansim.html). 
This process allows one to develop and test the observing strategy and analysis software in advance of instrument 
delivery. Furthermore it is possible to assess the impact of instrumental design aspects on the science objectives. 
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FIGURE 1. 7-element, SOOjum l-FX spaced hexagonal close-packed array geometry with 59" FWHM beams, and bolometer 
signals for a scan across a point source (path of source indicated by arrow). 
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The simulated sky, typically generated by convolving a synthetic catalogue of extra-galactic sources with the 
telescope beam [1], is assigned equatorial coordinates on the celestial sphere. The geometry of the array is defined 
and the beam positions on the sky are determined relative to the telescope bore-sight. The telescope is then arbitrarily 
located and the beams are scanned across the simulated sky model (in some pattern on an Alt-Az coordinate system 
for a given local time) and, using a rigorous calculation of the astrometry, the positional information and noiseless flux 
density time-series is determined for each bolometer in the array. 

More realistic time-series are produced by incorporating features characteristic of real instruments, the pointing 
performance of the telescope, and the non-negligible contribution of the atmosphere (sky noise and attenuation). 

Sources of instrumental noise in the time-series include: 1) addition of an uncorrected Gaussian system noise 
component; 2) multiplication by gain factors (or different responsivities) for each detector; 3) addition of independent 
low frequency l/f noise for each bolometer time-series to model the slowly-changing detector baselines in long 
duration scans; 4) convolution with an e _f / x impulse function to mimic the instantaneous response of the detectors. 
The instrumental noise components for a linear detector are summarized by the following expression for the response: 

^measured — G X F -\r I dark ( 1 ) 

where G is the detector gain, F is the incident flux and Idark is the dark current consisting of both broad-band and l/f 
noise components. The perfect astrometry may also be corrupted by adding random and systematic pointing errors. 
Fig. 1 shows the time-series for a 7-pixel feed-horn coupled array scanned across a point source at 500^/m. 



SKY NOISE 




FIGURE 2. A 360 x 180 degree 2 simulation of sub-mm sky emission at 850/jm from cells of water vapour distributed in a flat 
plane at an altitude of 600m above the telescope. 

Sky noise at sub-mm wavelengths is due to variable emission from cells of water vapour moving across the field-of- 
view of the telescope. This sky signal is significantly greater than the astronomical signal, and errors introduced 
through its imperfect subtraction lead to increased noise and artifacts in the reduced data. Thus, this important 
component of noise must be considered in all realistic simulations of ground-based observations. 

Data from the SCUBA camera [2] on the James Clerk Maxwell Telescope on Mauna Kea suggests the sky generates 
1 / /-like noise that is correlated across the array. An empirical model that reproduces the observed spectrum of signals 
measured with the SCUBA array at 850/jm assumes the sky consists of cells of water vapour, generated by convolving 
Gaussian noise with a symmetrical ( 1 //) 1 function, where / is the 2-D spatial frequency, distributed in a flat-plane 
at some fiducial altitude, moving at the wind-speed above the telescope aperture (Fig. 2). This planar sky-model is 
passed across the simulated beams to produce a strong, variable sky signal. 




FIGURE 3. Simulations of one-second point-and-stare total power observations with a fully-sampled (0.5-FA.) 40 x 40 pixel 
monolithic array at 850/jm with 14" FWHM beams. In the case of perfect flat-fielding (left panel), two significant sources with a 
S/N of 25 and 12a are clearly observed against a typical sky-background gradient. In contrast, the fainter of the sources blends into 
the background noise when the gains of the bolometers are only known to a precision of ~ 1% (right panel). 



Each 850/jm pixel in the SCUBA2 array under good weather conditions will be illuminated by a sky background 
power of 8pW with a photon noise in a Is integration of 4.3 x 10~ 5 pW. While the total power of the sky may be 
removed by subtracting the median signal of the off-source sky bolometers, strong gradients still remain and can 
dominate the astronomical signals. DC sky measurements with SCUBA are used to estimate these gradients by 
transforming the time-series into a spatial signal via the recorded wind speed, assuming the emission is from a fixed 
altitude above the telescope, and scaling it accordingly for the different background loading of the SCUBA2 pixels. 
The estimated median change in the sky level across an 8 arcmin field of view is 7 x lCT 4 pW (about an order-of- 
magnitude greater than the photon noise). These estimates of the expected gradients and DC power levels allow the 
proper scaling and bias of the planar-cloud model used to generate the sky component of bolometer signals in the 
simulation. 

Effective subtraction of the sky in the case of total power measurements depends critically upon the knowledge 
of detector gains and drifts. In Equation (1), G and Idark vary from pixel to pixel and need to be removed using 
separate operations. The gain G is slowing-changing, and characterized by performing flat-field observations with 
long integrations. The flat-fielding error must be smaller than the photon noise measured in the longest astronomical 
integration. In addition G is multiplied by the large sky signal (3> astronomical signal), small errors in its measurement 
will lead to residual artifacts upon dividing it out of the response. An initial estimate suggests that G will need to be 
measured to an accuracy of <C 1 %. Idark consists of broad-band (white) noise plus a slowly varying 1// component; the 
drift may be subtracted by measuring "dark" frames (Is integrations) in the absence of sky signal every few minutes. 

Fig. 3 demonstrates, for pointed total-power observations, two examples of the effects described above. In conclu- 
sion, the requirement to flat-field at sub-mm wavelengths with a precision of <C 1 % presents a major obstacle to the 
efficient use of large-format monolithic arrays on telescopes operated in total power mode. 
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